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ABSTRACT 

We analyzed ROSAT PSPC images of the bright, nearby {z < 0.1) galaxy clusters 
A401, A1656 (Coma), A2029, A2199, and A2256 to constrain their intrinsic shapes 
and baryon fractions; the intrinsic shapes of these clusters were analyzed previously 
by us using Einstein data (Buote & Canizares 1992; Buote 1992). Following Buote 
& Tsai we probed the aggregate structure of the clusters on scales ~ l.ShgQ Mpc to 
reduce effects of possible substructure on smaller scales a few hundred kpc). The 
ellipticities of the X-ray isophotes are typically e ~ 0.15 — 0.25 and display negative 
radial gradients highly significant at the 95% confidence limit for all the clusters except 
Coma. By assuming hydrostatic equilibrium and a variety of mass models we obtain 
ellipticities (mass ~ 0.40 — 0.55 for isothermal models of the total gravitating matter 
of the clusters; the (mass constraints change by < 10% upon consideration of the small 
temperature gradients shown by ASCA to be typical for rich clusters; the observed 
X-ray ellipticity gradients require that mass models with constant ellipticity have steep 
density profiles, p cx r""*. Estimates of the gas masses are highly insensitive to the 
ellipticities of the X-ray isophotes. The clusters in our sample have increasing fractions 
of gas mass to total mass with radius and have Mgas/Mfot = (4% — ll%)hgQ within 
a radius l.ShgQ Mpc, in excellent agreement with the results of White & Fabian and 
the Baryon Catastrophe proposed by White et al. (1993). Finally, the ellipticities of 
the dark matter distributions are essentially identical to (mass and are consistent with 
the shapes of dark halos predicted by N-body simulations and the shapes of the galaxy 
isopleths in the clusters in contrast to our previous conclusions using Einstein data. 
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1. Introduction 

The structures of galaxy clusters probe cosmological theories in a variety of ways. The 
influence of early tidal distortions in the formation of clusters (Binney & Silk 1979), the degree 
of dissipation of the dark matter (e.g., Strimple & Binney 1979; Aarseth & Binney 1978), and 
the density of the universe (Eisenstein & Loeb 1995) all affect the intrinsic shapes of clusters. 
Morphological differences between the X-ray-emitting gas and the underlying mass in clusters 
may reflect the nature of the dark matter itself (Kaiser 1991). Finally, the large baryon fractions 
observed for Coma (Briel, Henry, & Bohringer 1991) and other clusters (White & Fabian 1995) 
may have serious repercussions for many cherished notions of the standard cosmology (White et 
al. 1993). 

We (Buote & Canizares 1992 - hereafter BC; Buote 1992) have previously analyzed the shapes 
of a small sample of bright, low-redshift [z < 0.1) Abell clusters having no obvious subclustering 
(with the possible exception of Coma) using the Einstein Imaging Proportional Counter (IPC; 
Giacconi et al. 1979) and concluded that the ellipticities of the gravitating mass (e ~ 0.3) were 
smaller than the ellipticities of the galaxy isopleths (e ~ 0.5). However, the increasing evidence 
for substructure in clusters (e.g.. West 1995) calls into question any analysis of clusters that 
approximates them as relaxed systems. Buote & Tsai (1995a; hereafter BT) used the N-body / 
hydrodynamic simulation of Katz & White (1993) to test the reliability of X-ray constraints of 
intrinsic cluster shapes considering the effects of substructure. For the X-ray method to be reliable 
they concluded it is necessary that there is no subclustering on the same scale used to compute 
the shape of the cluster. Speciflcally, BT demonstrated that X-ray constraints of the aggregate 
shape of the Katz & White cluster on scales r ~ 1.5 Mpc are very insensitive to subclustering on 
scales <^ few hundred kpc. This is simply a statement that in a hierarchical clustering scenario the 
bulk of a cluster within r ~ 1 — 2 Mpc may be essentially relaxed even though a small subcluster 
(or two) has recently fallen in to the central regions. 

Following the suggestion of BT, in this paper we use data from the ROSAT satellite to probe 
the aggregate shapes of the bright, low-redshift [z < 0.1) Abell clusters A401, A1656 (Coma), 
A2029, A2199, and A2256 corresponding to clusters we analyzed previously with IPC data (BC: 
Buote 1992); some of these clusters have now been shown to have signiflcant subclustering in their 
cores (e.g., Briel et. al. 1991; White, Briel, k Henry 1993; Mohr, Fabricant, k Geller 1993). The 
image reduction and spatial analysis is described in §2. The hydrostatic modeling of the clusters, 
the results for intrinsic shapes, and the baryon fractions are discussed in §3. We present our 
conclusions in §4 

2. Spatial Analysis of X-ray Data 

The cluster images were obtained from the ROSAT Public Data Archive operated by the 
HEASARC-Legacy database at Goddard Space Flight Center; see Triimper (1983) for a description 
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of ROSAT and Aschenbach (1988) for a discussion of the X-ray telescope. We selected images 
observed with the Position Sensitive Proportional Counter (PSPC; Pfeffermann et al. 1987) 
instead of the High Resolution Imager (HRI) because of the PSPC's superior sensitivity, larger 
field of view, and, most importantly, all of the clusters were observed with the PSPC. Each cluster 
has at least one observation pointed near the center of the cluster, while a few of the clusters 
have multiple pointings, some of which are offset from the centers of the clusters. Only certain 
pointings, however, are useful for our analysis of the intrinsic shapes and mass distributions in the 
clusters. 

Interior to the 40' diameter ring of the window support structure of the PSPC the spatial 
resolution is substantially higher than outside (Hasinger et al. 1994). Outside the ring the PSPC 
"spokes" supporting the window fan radially outward from the ring and thus add structure to 
the PSPC images. Since the poor resolution and "spokes" outside the ring degrade and bias 
measurement of the ellipticity of the surface brightness (i.e. quadrupole; see §2.3) we restrict our 
analysis to regions interior to the PSPC ring (although see §2.2 regarding A2199). Moreover, 
when considering multiple pointings, we only used those that completely encircled a cluster about 
its center. In Table 1 we list the ROSAT sequence numbers and exposure times for the relevant 
observations of the clusters. 



2.1. Image Reduction 

To prepare the images for analysis we (1) removed time intervals of high background, (2) 
selected PI bins corresponding to photon energies between 0.5 and 2 keV, (3) corrected for 
exposure variations and telescopic vignetting, (4) merged multiple pointings for relevant clusters 
into one image, (5) subtracted the background (only for radial profile), and (6) rebinned the image 
into pixels corresponding to ~ 50/igQ^ kpc (Hq = SOhgo km s~^ Mpc~^). All of these reduction 
procedures were implemented with the standard IRAF-PROS software. 

All of the pointed observations were partitioned into many short exposures in order to 
maximize the observational efficiency of the ROSAT observing program. We examined the 
background light curves of the images for short-term enhancements indicative of contamination 
from scattered X-rays, especially from the sun, the bright earth, and the SAA. Only images for 
A1656 and A2256 required time-filtering. We list the effective exposure times of these filtered 
pointings in parentheses in Table 1. 

To minimize the effects of the X-ray background and the width of the PSPC point spread 
function (PSF; see §2.2) we selected photons only from energy channels between 0.5 and 2 keV. 
In addition, we rebinned the images into more manageable 15" pixels corresponding to 512 X 512 
fields. This pixel scale is the same as the exposure maps provided with the standard analysis 
systems software (SASS); note that the true resolution of the exposure maps actually corresponds 
to 30" pixels. 
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The images were then flattened using the SASS exposure maps. When dividing the images 
by these exposure maps, we corrected for both exposure variations and telescopic vignetting. In 
principle this correction depends on the energy of each individual photon, but for energies above 
0.2 keV the energy dependence is small and we neglect it (Snowden et al. 1994). 

We merged multiple pointings for A401, A1656, and A2256. Point sources common to the 
images of each cluster were used to align the flelds. A401, which has two images nominally 
centered on the same position, required a N-S shift of 1.8 it 0.5 pixels, slightly larger than the ~ 1 
pixel uncertainty (i.e. 15") expected of the pointings of ROSAT. The merged images for A1656 
included ~ SOOhg^ kpc within the ring of the PSPC. The flve pointings for A2256 easily flt 1.5 
hgQ Mpc radius within the PSPC ring. 

The next step is to remove point sources embedded in the cluster continuum emission. It 
is imperative to carefully remove such sources to prevent contamination of measurements of 
the ellipticity of the surface brightness (see BT). Our preferred method to remove sources is by 
"symmetric substitution" as outlined in BT. This method exploits the property that since our 
models used for analysis of the intrinsic shape (see §3.1) assume elliptical symmetry we may 
replace a localized region of the surface brightness with the corresponding regions related by 
reflection over the major and minor axes of the elliptical isophotes. Unfortunately, this method 
may only be implemented for an image containing very few sources and having well-deflned 
elliptical symmetry axes. Only A2029 and A2199 meet these requirements for the clusters in 
our sample. For the other clusters we removed sources by flrst selecting a source by eye and 
choosing an annulus around the source to estimate the local background. We then flt a second 
order polynomial surface to the background and replace the source with the background. This 
procedure is well suited for estimating quadrupole moments of high S/N cluster images; see Buote 
& Tsai (1995c) for a thorough discussion. 

We estimated the background for each cluster from source-free regions outside the ring of the 
PSPC flelds; these regions were away from the held center at radii r ~ 40' for A401, 30' for A2029, 
45' for A2199, and 55 — 60' of the cluster center for A2256. The background rates in the 0.5 - 2 
keV band of these regions are listed in Table 1. For Coma we adopted the value of White, Briel, 
& Henry (1993) since they explored regions much farther from the cluster center than we do. We 
mention that our background estimate for A2256 agrees well with the value 2.24 X 10""* cts s~^ 
arcmin"^ of Briel & Henry (1994). 

Finally, we rebinned the images of each cluster in order to probe the aggregate cluster 
structure as outlined in BT. We rebinned the images into pixels corresponding to ~ 50/igQ^ kpc for 
computation of the radial proflles of the clusters; this arbitrary choice provides enough resolution 
to determine the radial structure suitable for analysis of the aggregate radial structure of the 
cluster. However, because reliable computation of the ellipticity (see §2.3) requires a sufficiently 
large number of pixels in a given aperture, we generally selected a finer pixel scale for computing 
the ellipticities (see Table 2). 
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Contour plots of the reduced images for each cluster are plotted in Figure 1 using the pixel 
scales of the radial profiles listed in Table 2. The image for Coma stands out since we are viewing 
a much smaller region (r ~ SOOhg^ kpc) than for the other clusters (r ~ ISOOhg^ kpc). Binned 
into these large pixels Coma appears very smooth whereas on smaller scales there is significant 
structure (see §2.2). This is in keeping with our strategy to measure the aggregate structure of the 
clusters to reduce possible nonequilibrium effects of small-scale subclustering. 

2.2. Radial Profile 

We constructed the azimuthally averaged radial profile for each cluster using the background- 
subtracted images. First we computed the centroid of the cluster within a circular aperture 
containing ~ 90% of the total fiux. This centroid was determined by selecting an origin by eye 
and then iterating until the centroid changed by < 0.1 pixels. For the clusters with known cooling 
fiows (i.e. A2029 and A2199) we excluded the region interior to the radius where the cooling time 
~ cluster age; this radius is typically ~ 100 kpc (see Sarazin, O'Connell, & McNamara [1992] 
concerning A2029; see Fabian [1994] for a review of cooling fiows in clusters). We exclude the 
cooling fiow region since our models in §3.1 do not apply there. 

When computing the centroid of A2256 we must carefully consider the well-known substructure 
in the core (Briel et al. 1991). Our analysis of the mass distribution in §3.1 assumes the cluster is 
in hydrostatic equilibrium. We prefer to exclude the region in A2256 were the gas is unrelaxed. 
To determine where in the cluster the gas is likely to be relaxed we appeal to the power ratios of 
A2256 (Buote & Tsai 1995b, c). The power ratios classify clusters according to their dynamical 
state and depend on the scale of the cluster being probed. Buote & Tsai (1995c) find that on 
a scale of O.Shg^ Mpc A2256 has power ratios similar to clusters that are manifestly unrelaxed. 
However, on LOhg^ Mpc scales A2256 has power ratios similar to clusters with only a small 
amount of substructure and is morphologically much closer to evolved clusters like A2029 than 
on the O.ShgQ Mpc scale. This effect of more complete relaxation away from the subcluster is 
supported by the temperature distribution (Briel & Henry 1994) which is complex (with measured 
azimuthal variations) within O.Shg^ Mpc and nearly uniform outside l.Ohg^ Mpc. To achieve a 
balance between eliminating as much of the inner ~ O.ShgQ Mpc because of substructure and 
sampling enough of the radial profile to usefully constrain the mass models in §3.1, we decided 
to excise the inner ~ 300hgQ kpc from analysis which effectively encompasses the subcluster in 
question. We thus compute the centroid in an annular aperture from (300 — 1500)/igQ^ kpc. 

For all of the clusters we computed the radial profiles out to ~ l.ShgQ Mpc where the specific 
outer bin is determined by a S/N criterion (see below). Both Coma and A2199, however, fall well 
short of l.ShgQ Mpc within the PSPC ring. We decided to extend the radial profile of A2199 
outside the ring in order to compute its baryon fraction to larger distances (see §3.3); the radial 
profile is not as sensitive as the ellipticity to the larger PSF and spokes outside the ring. Extending 
the radial profile does not improve our constraints on the intrinsic shape of A2199 since it is the 
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X-ray ellipticity that is the most important determinant which is still confined to within the ring 
(see Buote & Canizares 1995a). We do not extend the radial profile of Coma since its baryon 
fraction has been analyzed in detail by others (White et al. 1993). 

We rebinned the radial profiles so that all bins had an appropriate minimum signal-to-noise 
ratio {S/N): 7 for A2029 and A2199, and 20 for A2256; A401 and Coma required no additional 
binning. The background-subtracted radial profiles are displayed in Figure 2. All of the radial 
profiles appear to be quite smooth for our chosen bin sizes. In particular, Coma's profile displays 
no large irregularities when binned into ~ 50hgQ kpc pixels, but is known to be very lumpy when 
probed on galactic scales of ~ 5 kpc (e.g.. White et al. 1993; Davis & Mushotzky 1993). This 
property that the cluster is sufficiently relaxed for analysis of its aggregate structure while it is 
clearly not a smooth, equilibrium configuration on smaller scales exemplifies the arguments of BT 
and is the driving force for our present investigation. 

A convenient parametrization of the X-ray radial profiles of galaxy clusters is the (3 model 
(Cavaliere & Fusco-Femiano 1976; Jones & Forman 1984; Sarazin 1986), 



where and (3 are free parameters. The (3 model is especially useful as (1) a benchmark for 
comparison ofT,x{R) to other clusters and (2) an analytic parametrization of the radial parameters 
for computing the gas mass. In order to obtain physical constraints on the parameters ax and (3, 
we convolved Ex with the off-axis PSPC PSF (Hasinger et al. 1994) and performed a fit to the 
radial profile; note for the evaluation of the PSF we set the energy to 1 keV and for A1656 and 
A2256, which have pointings on different regions of the clusters, we fixed the off-axis angle to 10'. 
Since our radial bins are larger, in some instances substantially larger, than the width of the PSF, 
including the PSF in the fits has a negligible (^ 1%) effect on the fitted parameters. Nevertheless, 
we include the PSF in the fits for completeness. 

In Table 3 we list the best-fit parameters and 95% confidence limits on two interesting 
parameters; the best-fit models are also plotted in Figure 2. The (3 model is an excellent visual 
fit to the radial profiles for all of the clusters, although the values are rather large. For the 
purposes of analyzing the aggregate structure of the clusters on ~ 1 Mpc scales, though, the fits 
give very acceptable descriptions of the surface brightness distribution. Except for A2256, the 
fitted parameters are generally insensitive to the aperture size used to evaluate the centroid; i.e. 
'^(ix ^ 5%, A/3 ^ 1%. A2256, which is known to have a large centroid shift (Fabricant, Rybicki, 
& Gorenstein 1984; Mohr et al. 1993), has substantially different parameters depending on the 
chosen aperture used for computing the centroid. The values listed in Table 3 for A2256 are for the 
centroid evaluated in an annular aperture (300 — 1100) hg^kpc representing ~ 90% of the fiux in 
the annulus (300- 1500) /igo^kpc. If we instead use a circular aperture of radius ISOOhg^ kpc then 
we obtain best-fit parameters = 5'. 21 and (3 = 0.815. These values are in excellent agreement 
with those obtained by Briel & Henry (1994) who removed the effects of the subcluster by excising 




2-1 -3/3+1/2 



(1) 



- 7 - 



a pie slice from the image with origin located at the primary cluster center and edges extending 
throughout the whole cluster. Our fitted /3-model parameters when including the subcluster agree 
with those of Briel & Henry when they exclude the subcluster in a pie slice; i.e. excluding the 
subcluster itself has a minimal effect on the radial profile but when the whole interior is excluded 
the /3-model parameters are noticeably affected. This may simply refiect the steepening of the 
radial profile outside the core. As stated above, we prefer to exclude from analysis the region 
r < 300/i8o^ kpc because the substructure there may signal substantial departures from hydrostatic 
equilibrium. Hence, except for A2256, the clusters have and (3 that agree well with those found 
in the literature (e.g., Jones & Forman 1984; Briel, Henry, & Bohringer 1991). 

2.3. EUipticity 

We quantify the elongation of the clusters using the iterative moment technique introduced 
by Carter & Metcalfe (1980). This method is particularly suited to measure the aggregate 
ellipticity of the clusters in a large aperture (see BT and Buote & Canizares [1994] for detailed 
discussions). In essence, this technique entails computing the analogue of the two-dimensional 
moments of inertia arrived at by iterating an initially circular region; the square root of the ratio 
of the principal moments is the axial ratio and the orientation of the principal moments yields 
the position angle. The parameters obtained from this method, and 9m, are good estimates 
of the ellipticity (e) and position angle {9) of an intrinsic elliptical distribution of constant shape 
and orientation. For a more complex intrinsic shape distribution, and 9m are average values 
weighted heavily by the outer parts of the region. 

We applied a simple Monte Carlo procedure to characterize the uncertainties on ^m and 9m 
due to undetected point sources and Poisson noise. For each cluster we simulated (3 models having 
the best-fit and (3 obtained in §2.2 modified to have a constant ellipticity and orientation; 
i.e. we replaced R = ^/x^ -\- with the elliptical radius \Jx^ ^ y"^ /q^^ where q is the constant 
axial ratio. These models were scaled to have the same number of counts (or count rate) as the 
background-subtracted PSPC images of the clusters. We then added the uniform background 
for each cluster (see Table 1). To each image we added point sources having spatial properties 
consistent with the PSPC PSF and numbers consistent with the logA^(> S) - \ogS distribution 
given by Hasinger (1992); see Buote & Tsai (1995b) for a detailed discussion. Finally, we added 
Poisson noise to the composite images. Since tM is unaffected by a uniform background (see 
Carter & Metcalfe) we did not then subtract the background. We performed 1000 simulations 
each for a suite of input ellipticities (e^; = 1 — q) for each cluster. 

To determine the 95% confidence intervals on the measured tM we proceeded as follows; 
for now we focus our attention on a particular aperture size. We arrange the results of the 
1000 simulations for a given input t^ into ascending order of measured ellipticity tM] i-e. 

^ < ■ ■ ■ < fM°°- '^^^ 95% upper limit for this model is defined to be the value of tM 
corresponding to the 0.95 X 1000 = 950th value of tM in the ordered array; i.e. ej||°. The 95% 
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confidence lower limit of the real data is given by the model with input whose e]|| just equals 
the measured value of from the real image (see below). Similarly, the 95% confidence upper 
limit of the real data is given by the model with input whose just equals the measured value 
of eM from the real image. 

We estimated the confidence limits for the position angles 9m from the model having = 
of the data for the specific aperture size. We define the 95% confidence limits to be (6j^,6^^). 
This procedure is not the most rigorous means to determine the uncertainties on 9m because 
the value of ^m in each simulation is not equal to e^- Since, however, we do not use 9m in our 
modeling (see §3.1) the estimates are sufficient for our purposes. 

In Table 4 we list the values of ^m and 9m and their associated 95% confidence limits for 
different apertures (in 5 pixel increments) for each of the clusters. Similar to our analysis of the 
radial profiles, we excluded the cooling-fiow regions of A2029 and A2199 from the ^m computation. 
Similarly, we excluded the region a ^ 300hgQ kpc for A2256 to lessen the infiuence of the core 
substructure. Generally the values of ^m are larger than the ellipticities obtained by BC (Buote 
1992) from Einstein IPC images of these clusters; note BC used a slightly modified version of 
the iterative moment technique. Other published ellipticities for these clusters using IPC data 
are consistent with the results of BC (Buote 1992). For example, using different techniques to 
measure the ellipticity Fabricant et al. (1984), Buote (1992), Mohr et al. (1993), and Davis & 
Mushotzky (1993) measured an ellipticity of ~ 0.2 for A2256 at 10' - all significantly less than the 
value of ~ 0.3 computed in this paper. For Coma BC and Mohr et al. obtained ellipticities at 
15' in good agreement with the value of 0.16 obtained in this paper although Davis & Mushotzky 
obtained ~ 0.27. The good agreement between IPC and PSPC results for Coma is expected since 
the effect of the PSF should be least important for Coma of all the clusters in our sample. Hence, 
the discrepancy between the ellipticities of the clusters computed with IPC data and the results 
of the present paper refiects the superior spatial resolution of the PSPC. 

All of the clusters except A1656 have ^m that are consistent within their uncertainties with 
being monotonically decreasing with increasing aperture size; in fact, the ellipticities of Coma 
appear to be increasing with radius but over a much smaller range. The position angles vary 
^ 10° over the regions probed except for Coma and A2029 which have ~ 15° variations. The 9m 
are generally in good agreement with the previous IPC studies and with the orientations of the 
galaxy isopleths (Carter & Metcalf). 

Following BT we designated two apertures within the ~ LShg^ Mpc regions for use in 
constraining the mass models in §3.1; for Coma we only used the largest aperture available. 
Using two large, well-separated apertures is in keeping with our scheme to analyze the aggregate 
structure of clusters while still allowing us to obtain information on any ellipticity gradients. We 
denote these special apertures with an asterisk in Table 4. 



- 9 - 



3. Intrinsic Shapes and Mass Profiles from the X-rays 

3.1. Method 

The technique we use to constrain the intrinsic shape and mass profile of the total gravitating 
mass for a cluster of galaxies is discussed in detail by BT. Buote & Canizares (1994,1995a) describe 
how we use the the observed light distributions to infer the shape and profile of the dark matter 
from X-ray images. We refer the reader to these papers for exposition of the modeling procedures 
we employ in this paper. Here we summarize only a few small additions to the method introduced 
for this particular study. 

Because of the increasing evidence from gravitational-lens observations of giant arcs (e.g., 
Miralda-Escude & Babul 1995) and N-body simulations (e.g., Navarro, Frenk, & White 1995) that 
galaxy clusters have small or perhaps no core radii, in this paper we will explore density profiles 
that are simple generalizations of the Hernquist (1990) model, 

Ph cc (am)""^ (i?c + am)~^ , (2) 

where Rc is a core parameter, a is the semi-major axis of the bounding spheroid of the SMD 
(Spheroidal Mass Distribution - see BT), and m is the dimensionless spheroidal radius. For 
am <^ Rc, ph ~ and thus these densities have no core. At large radii, am ^ Rc, pu ~ r"^""'"^) 
and falls as r~'^ for the Hernquist (1990) density (ra = 3) we used in previous studies. We consider 
models having n ranging from 1 — 3 which covers plausible ranges of density profiles of clusters. 
These models complement the "cored" power-law SMD models [p oc [Rc + {am)'^]~^) that we also 
use to model the masses of clusters. 

To describe the mass distribution of the X-ray-emitting gas we use simple models having 
isodensity surfaces that change shape with radius. We are lead to this choice since the X-ray 
ellipticity profiles of the clusters measured from the PSPC data (Table 4) display significant 
gradients. An accurate model for the gas distribution is desirable not only for obtaining a precise 
measure of the total gas mass, but also for constraining the shape and profile of the dark matter. 
That is, in galaxy clusters the dominant mass component is the dark matter which typically 
contains 80% — 90% of the mass. The gas mass is 10% — 20% of the total while the mass of the 
stellar material, which is not nearly as well constrained as the gas, is generally much smaller than 
the gas mass (e.g., the case of A2256 - Briel & Henry 1994). Since the total mass of the galaxies 
is not well constrained, and in any case probably insignificant with respect to the gas, we neglect 
it (i.e. Mfot = Mgas + M^m)', note that for elliptical galaxies the opposite is true - the gas may 
be neglected in favor of the stellar mass (e.g., Buote & Canizares 1994,1995a). 

The models we use to describe the gas mass are the oblate (and prolate, see below) spheroids 
of varying ellipticity introduced by Ryden (1990). Ryden computes the monopole and quadrupole 
terms of the gravitational potential generated by a mass distribution of concentric, oblate, 
spheroidal shells whose eccentricity (e^ = 1 — q^) varies as a function of the semi-major axis of the 
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shell. Like Ryden, we consider spheroids that have eccentricity, 

^^^^ = TTW' 

where we have inserted the extra parameter cq to allow for an arbitrary central eccentricity. 
This simple function yields surprisingly good descriptions of the ellipticity profiles of the X-ray 
isophotes (see Table §4). Since we also desire to explore prolate models we extend Ryden's 
formalism to the prolate case in the Appendix. 

For each cluster we construct Ryden models for the gas mass. Since (3 models describe the 
radial profiles of the clusters in our sample so well (§2.2) we model the three-dimensional density 
profiles as pgas = (a^ + a^)""^^/^. Here and f3 are the best-fit parameters obtained from fitting 
the (3 model convolved with the PSPC PSF to all of the radial bins for each cluster. To completely 
specify the gas model we assign e{a) (eqn. [3]) to pgas appropriate to the ellipticity of the X-rays. 
That is, we generate a model X-ray surface brightness by projecting p'^^^ onto the sky plane and 
convolving it with the PSPC PSF. By comparing the ellipticities of this model with those of the 
real image (see Table 4) we determine the parameters of e{a). 



3.2. Shapes of the Composite Mass Distributions 

Our procedure to constrain the hydrostatic models using the radial profiles (§2.2) and 
ellipticities (asterisked values in Table 4) is described in detail in BT §4.2. For all of the spheroidal 
models we set the symmetry axis to lie in the plane of the sky; i.e. we are not attempting to 
uncover projection effects in this analysis, although we believe them to be small for these clusters 
(see §3.4). We set the semi-major axis of the SMD models to l.ShgQ Mpc for all the clusters. 



3.2.1. Isothermal Models 

Recent ASCA observations show that that gas is nearly isothermal within the central ~ 1 
Mpc of rich clusters (Mushotzky 1995). Since the shape of a cluster inferred from X-ray analysis 
is mostly insensitive to small temperature gradients (BT; also Strimple & Binney 1979 and Buote 
& Canizares 1994), we focus our attention on isothermal models. The effects of temperature 
gradients consistent with PSPC data are discussed in §3.2.2. 

We list in Table 5 the 95% confidence limits on the ellipticities of the composite mass {tmass) 
for the isothermal models that best characterize the X-ray data. Similar to our experience from 
fitting the [i models (§2.2), the isothermal models give excellent visual fits to the data but with 
rather large values of ■ The two cooling-fiow clusters, A2029 and A2199, require steep mass 
density profiles {pmass ~ r~'^) to reproduce the observed ellipticity gradients of the X-rays; even 
with this steep density the models for A2029 only marginally reproduce the X-ray ellipticity 
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gradient within the estimated 95% confidence uncertainties. The ellipticities in Table 5 for A2029 
and A2199 are for the Hernquist (1990) models (i.e. ra = 3 in eq.[2]), although the cored power-law 
models with r] = 2 give essentially identical values. 

The clusters A401 and A2256, which have no observed cooling fiows (Fabian 1994), also 
have steep X-ray ellipticity gradients that imply steep mass density profiles {pmass ~ r""*). For 
these clusters, though, we list the results in Table 5 for the cored power-law models with r] = 2 
since their radial profiles have large core radii (§2.2); as was the case for A2029 and A2199, the 
ellipticities for A401 and A2256 are not very sensitive to the density model and we obtain similar 
results for the Hernquist models. Coma, for which we only use a single aperture to constrain 
the ellipticity (and thus have no information on the large-scale ellipticity gradient), is fit equally 
well by all of the models considered. We show in Table 5 the results for both the r] = 1,2 cored 
power-law models {pmass ~ r~'^ and p^ass ~ r~*). 

The ellipticities required by the isothermal models are quite large for all of the clusters, much 
larger than those we obtained from analysis of Einstein IPC data in BC. However, the ellipticities 
for A2256 are in good agreement with IPC analyses of Fabricant, et al. (1984) and Buote (1992). 
We attribute the discrepancy to BC not accounting for the steep ellipticity gradients of the X-ray 
isophotes obtained for these clusters (see Table 4). In BC we selected only one aperture size to 
constrain the X-ray ellipticity and the size of the aperture was only a few times the width of the 
IPC PSF for all the clusters but Coma (and Perseus). Since in BC we did not incorporate the 
IPC PSF into the hydrostatic modeling procedure we did not accurately estimate the uncertainty 
due to possible ellipticity gradients being smeared out by the IPC PSF. In the case of Coma the 
discrepancy is due primarily to the larger aperture we use in this paper to compute the X-ray 
ellipticity; when we use the same aperture size as in BC we obtain (mass in good agreement. The 
agreement of the shape of A2256 with previous analyses appears to have been due to a fortuitous 
conspiracy of the different region analyzed in this paper and the larger X-ray ellipticities due to 
the better spatial resolution of the PSPC. 

3.2.2. Models with Temperature Gradients 

We now consider temperature gradients allowed by the PSPC data. Unfortunately, since the 
PSPC energy band corresponds to much lower energies than the temperatures of the rich clusters 
in our sample we are unable to obtain precise constraints on the temperature gradients. In general 
we find that the uncertainty in the gradients allowed by the PSPC data dwarf those obtained 
by PV-phase ASCA observations of rich clusters (Mushotzky 1995). Instead of giving serious 
consideration to models that do not appear to be realistic, we will simply illustrate this uncertainty 
with the case of A2029. We defer serious consideration of the effects of real temperature gradients 
in these clusters on the shapes to analysis of ASCA data (Buote & Canizares 1995b). 

In order to obtain the narrowest constraints on the temperature gradient we placed plausible 
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restrictions on the spectrum and the temperature profile. Our intention is to demonstrate that 
even when such restrictions are applied there is still an implausibly large range of temperature 
gradients allowed. We first divide the surface brightness of A2029 into two annular regions 
having equal S/N: (140-336,336-1400) hg^ kpc; we exclude the region r < lAOhg^ kpc because 
of the cooling fiow. The temperature of each region was determined by fitting (with XSPEC) a 
Raymond-Smith (1977) model with Galactic absorption and 50% solar metallicities. We obtain 
= 7.7(6.7 - 9.3) keV and Tout = 6.1(4.8 - 8.1) keV at 68% confidence. 

We take the three-dimensional temperature profile to follow a broken power law, 

T{r) = To r<ro (4) 

T{r) = To(^-^y r>ro (5) 

where ro = 140/igQ^ kpc, and Tq and p are free parameters. For the purposes of the determination 
of the temperature gradients implied by the spectral constraints, we suppress the information 
supplied by the elongation of the X-ray isophotes which only has a small effect on the parameters 
To and p. However, for the hydrostatic models (see below) we force the isotemperature surfaces 
given by equation [5] to be stratified on the isopotential surfaces. 

We construct the emission-weighted temperature {{T)) projected along the line of sight (dl) 
into an area dA, 

= ' 

where jgas is the gas volume emissivity and T,^ = /f^ jgas dl is the surface brightness. To 
determine Tq and p we evaluate equation [6] over A corresponding to the two annular regions 
(140-336,336-1400) hg^ kpc and then perform a to the 68% confidence limits on the 

temperatures; note we defined the temperatures in the bins to be the middle of the 68% confidence 
temperature results, and the la weights are half the widths of the confidence intervals. 

The best-fit parameters are Tq = 8.8 keV and p = —0.27; the 68% lower limit is Tq = 11.1 
keV and p = —0.68 and the upper limit is Tq = 6.7 keV and p = 0.07. The slopes of the best-fit 
(p = —0.27) and 68% lower limit {p = —0.65) temperatures imply substantially larger gradients 
than observed in real clusters (Mushotzky 1995), although the 68% upper limit slope is modest 
{p = -|-0.07). We list in Table 6 the derived mass ellipticities for the temperature gradients 
determined above using the same mass density models for A2029 given in Table 5. As expected, 
the ellipticities for the small gradient given by the 68% upper limit are shifted down only by 
~ 0.03 with respect to the isothermal results; even (mass for the best-fit temperature parameters 
are shifted upwards by only ~ 0.05. However, the very large gradient given by the 68% lower limit 
temperature gradient implies implausibly fiat mass distributions for a non-rotating, self-gravitating 
spheroid (Merritt k Stiavelli 1990; Merritt k Hernquist 1991). 

In the previous discussion we have only considered the effects of radial temperature gradients 
on the derived shape of the total gravitating cluster mass. Azimuthal temperature variations 
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can distort the shapes of the X-ray isophotes since the X-ray emissivity also depends on the gas 
temperature via the plasma emissivity (convolved with the spectral response of the detector). 
Since the PSPC energy pass band corresponds to smaller energies than the typical temperatures 
of rich clusters the plasma emissivity convolved with the spectral response of the PSPC is 
effectively constant for the clusters in our sample (see NRA 91-OSSA-3, Appendix F, ROSAT 
mission description); i.e. the isophote shapes measured with the PSPC are largely insensitive to 
temperature variations. Of more serious concern is that large azimuthal temperature gradients 
likely imply strong departures from hydrostatic equilibrium in which case a non-equilibrium 
analysis is required; e.g., the temperature distribution of the core of A2256 has been modeled 
using hydrodynamic simulations (Burns 1995). However, such clusters often exhibit substructure 
which allows one to judge (without knowing the temperature) where in the cluster hydrostatic 
analysis is appropriate (as in A2256 - see §2.2). 



3.3. The Total Gravitating Mass and Baryon Fraction 

In order to determine the shapes of the dark matter we must first know the ratio of the gas 
mass to total gravitating mass (see §3.1) for the clusters. To compute the gas mass we simply 
integrated the models for pg^s in §3.1. We normalized pg^s for each cluster by first constructing 
the X-ray emissivities, j'^asj (see eq.[9] of Buote & Canizares 1994) from pg^s and then projected 
jgas onto the sky plane to obtain Ti^- We then normalized Ti^ to the total fiux between 0.5 - 2.0 
keV determined by fitting the spectrum (with PROS) to a Raymond-Smith (1977, updated to 
1992 version) model with Galactic absorption (David et al. 1993 who use the results of Stark et al. 
1992) and temperatures from Edge et al. (1990); see Table 2. By normalizing to the fiux we have 
completely specified pg^s and hence the gas mass. Our estimates for the gas mass take into account 
the following uncertainties: (1) 95% confidence statistical errors on the /3-model parameters 
(see Table §3) used to define the radial profile of Pgas^ (2) oblate and prolate geometry, and (3) 
variation in fiux and plasma emissivity due to the 90% uncertainties in the gas temperatures from 
Edge et al.. There is an additional source of uncertainty due to the restriction of the adopted 
models for pg^s to having radial profiles of the [i models. Since, however, the [i models describe the 
clusters so well over the regions considered the uncertainty in the gas mass should be quite small, 
certainly ^ 10%. We list in Table 7 the gas masses computed in spheres of radii (0.5, 1.0, l.b)h'^Q 
Mpc considering (1) - (3) for all of the clusters; note for Coma we only compute the gas mass out 
to 0.5/i8o^ Mpc. 

The gas masses are extremely well constrained for the clusters, generally to better than a few 
percent. Our values for Mg^s appear to be slightly larger than estimates found in the literature. 
White & Fabian (1995) used Einstein IPC data to compute gas masses and total masses of clusters 
assuming spherical symmetry for a sample of Abell clusters, including A401 and A2029. They 
obtain Mg^s = 1.32 ± 0.07 X for A401 (r < 1.265 Mpc) and Mgas = 1.26 ± 0.11 X 

for A2029 (r < 1.291 Mpc), all quantities evaluated for Hq = 50 km s~^ Mpc~^; note that 
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Mgas cc Using the same apertures and Hubble constant as White & Fabian we obtain 

Mgas = (1-99 - 2.02) X 10^'^Mq for A401 and Mgas = (1-90 - 1.97) X 10^'^Mq for A2029. Although 
our values are only ~ 50% larger, the discrepancy appears to be significant within the estimated 
uncertainties. White & Fabian, however, acknowledge that their estimates of the gas masses are 
likely to be somewhat conservative because of the procedure they employ to obtain core radii 
for the gas profiles. If we neglect the ellipticity of the surface brightness the masses are hardly 
affected; e.g., Mgas = (1-85 - 1.92) X lO^"* Mq for A2029. This confirms earlier statements that 
ellipticity has negligible impact on estimates of gas masses in clusters (e.g.. White et al. 1994). 
By extrapolating our model for pgas to a 3 Mpc {Hq = 50 km s~^ Mpc~^) radius for A2029 we 
obtain Mgas = (5.6 — 5.9) X 10^"* Mq, in good agreement with Jones & Forman's (1984) value of 
Mgas = 5.3 X 10^"* Mq obtained from IPC data considering their expected (unstated) uncertainties. 

The gas mass for A2256 has been computed by Briel & Henry (1994) using the ROSAT PSPC 
data. They obtain Mgas = 6.33±1.17x lO^'^M© within a 1.4 Mpc radius {Hq = 100 km s"! Mpc"!). 
At this distance we obtain (neglecting the ellipticity of the gas) Mgas = (7.46 — 7.84) X 10^"* Mq, 
only ~ 15% larger than the result of Briel & Henry. This small discrepancy can likely be attributed 
to that fact that Briel & Henry exclude a large section of the data from analysis different from us 
(see §2.2) in order to avoid the subclump in the interior. 

Using the isothermal spheroidal models of Table 5 we present in Table 7 the total gravitating 
masses of the clusters computed in spheres of radii (0.5, 1.0, 1.5)hgQ Mpc. The uncertainties in 
the masses refiect 95% confidence statistical errors on Rc and n (or r]) from the mass models and 
the 90% errors on the temperatures from Edge et al. which is the dominant source of uncertainty. 
It is worth mentioning that because our mass models have steep density profiles (i.e. p ~ r""*) 
the results of Table 7 differ from the total masses obtained assuming the cluster is a /3 model; i.e. 
a (3 model corresponds to a logarithmic potential (e.g.. Appendix C of Trinchieri, Fabbiano, & 
Canizares 1986) and thus p ~ at large radii. The difference between models is highlighted by 
the case of A2199 for which Mushotzky (1995) obtains Mtot = (1-5 - 1.7)h^Q X 10^"* Mq within 
O.ShgQ Mpc using the /3 model - we obtain the same result using the derived /3-model parameters 
in Table 3. However, using the p ~ r""* models we obtain (in Table 7) a mass that is over a factor 
of six(!) larger than the /3 model at O.Shg^ Mpc, although the descrepancy between the models 
decreases with increasing radius. However, the discrepancy between models is not so pronounced 
for all the clusters, especially at large radii. When scaled to 1.4 Mpc {Hq = 100 km s~^ Mpc~^), 
we obtain Mfot = (12 — 21) X 10^"* Mq for A2256 using the p ~ r""* models, only slighly larger than 
the values we obtain assuming a /3 model Mfot = (7.0 — 10.8) X 10^"* Mq. The /3-model values also 
happen to be in excellent agreement with the result of Briel & Henry (1994), ~ 9.5 X 10^"* Mq, 
who also used the isothermal /3-model to compute the mass. 

As expected, we find that the masses derived assuming spherical symmetry agree very 
well with those incorporating the ellipticity of the X-ray isophotes; the viability of X-ray mass 
estimates assuming spherical symmetry has already been addressed in the literature (Tsai, Katz, 
& Bertschinger 1994; Navarro et al. 1995). All of the clusters in our sample that are measured to 
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large radii have Mgas/Mfot increasing with radius. Scaling to Hq = 50 km s ^ Mpc ^ the mass 
ratios are (8 — 22)% consistent with the results of White & Fabian (1995) for their sample of 19 
Abell clusters; Mgas/Mtot cc h'^^^. 

3.4. Shapes of the Dark Matter Distributions 

Using the ratios Mgas/Mfot in Table 7 we determine the ellipticity of the dark matter 
distributions, eoM, for the clusters. In order to consider the maximum possible effects of the 
self-gravitation of the gas we also consider the upper limits of Mgas/Mfot scaled to Hq = 50 km 
s~^ Mpc~^; i.e. this scaling amounts to doubling the values of Mgas/Mfot in Table 7. The results 
for tDM corresponding to the mass models in §3.2.1 (see Table 5) are listed in Table 8. The 
ellipticities of the dark matter are hardly affected upon consideration of the self-gravitation of the 
gas; i.e. for all the clusters toM ~ ^tot + 0.02. 

Also listed in Table 8 are the la ellipticities of the galaxy isopleths, tg^u computed by Carter 
& Metcalfe (1980) using the same iterative moment technique of §2.3. All of clusters have toM 
consistent with tg^i within the (1 — 2)a uncertainties of tg^u although toM for A2256 is only 
marginally consistent with the 2a lower limit for tg^i- The case for A2256 is actually uncertain 
since Fabricant, Kent, & Kurtz (1989) obtained a much smaller lower limit for tg^i = 0.23 (90% 
confidence) using a different sample with more redshifts and a different technique to compute the 
ellipticity from Carter & Metcalf. Hence, the shapes of the dark matter distributions and galaxy 
isopleths are consistent for each of the clusters in our sample contrary to our previous conclusions 
in BC. (We mention that tg^i for A1656, A2029, and A2199 obtained by Plionis, Barrow, & Frenk 
[1991] using a different technique also agrees with toM- Since Plionis et al. find that the intrinsic 
tgai distribution for clusters peaks at ~ 0.5, the clusters in our sample are likely viewed nearly 
edge-on and thus our conclusions for the dark-matter shapes may be typical.) 

Finally, the ellipticities of the dark matter distributions for the clusters are consistent with 
the predictions from dissipationless formation of halos in a universe filled with Cold Dark Matter 
(Frenk et al. 1988; Efstathiou et al. 1988). The dark matter halos produced by these N-body 
simulations are triaxial and generally have toM < 0.5 but several have toM approaching 0.67. Our 
mass models fit the observed X-ray radial profiles and ellipticities equally well for both oblate and 
prolate models which is consistent with the halos being typically triaxial. These results are also 
consistent with the cluster formed in the simulation of Katz & White (1993) which has toM = 0.54 
within l.bh'^Q Mpc (see BT). 

4. Conclusions 

We have analyzed ROSAT PSPC X-ray images of five bright, low-redshift [z < 0.1) Abell 
clusters for the purpose of constraining their intrinsic shapes and mass profiles. The intrinsic 
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shapes of the clusters in our sample (A401,A1656,A2029,A2199,A2256) were analyzed previously 
by us (Buote & Canizares 1992; Buote 1992) using Einstein IPC X-ray images. In this paper we 
specifically follow Buote & Tsai's (1995a) procedure to constrain the aggregate shapes of clusters 
on large scales (~ 1.5 Mpc) to reduce the effects of possible substructure on small scales few 
hundred kpc). 

We computed the azimuthally averaged radial profiles and ellipticities (i.e. quadrupole 
moments) of the X-ray surface brightness distributions within ~ l.ShgQ Mpc for the all clusters 
except Coma where we were restricted to ~ O.Shg^ Mpc. Fitting /3 models to the radial profiles 
yields very precise constraints on the core radii and (3 parameters in good agreement with previous 
results from Einstein (e.g., Jones & Forman 1984). The ellipticities of the X-ray images are also 
tightly constrained in their 95% confidence intervals and are systematically larger (e ~ 0.15 — 0.25) 
than those (e ~ 0.10 — 0.20) obtained by Buote & Canizares (1992; Buote 1992) using Einstein 
data. This discrepancy is a direct result of the superior spatial resolution of the PSPC. All of the 
clusters (except Coma) display substantial ellipticity gradients within r ~ LShg^ Mpc not seen in 
the Einstein data. 

Using the X-ray radial profiles and ellipticities and assuming hydrostatic equilibrium we 
constrained the ellipticity of the total gravitating matter following Buote & Tsai (1995a). 
Isothermal mass models for the clusters yield ellipticities, e^ass ~ 0.40 — 0.55 (95% confidence), 
which are systematically larger (i.e. Ae^ass ~ 0.15) than obtained by Buote & Canizares (1992; 
Buote 1992) using Einstein data; we understand this discrepancy to be a combination of effects 
due to different cluster regions being probed, modeling differences between the two investigations, 
and most importantly the X-ray ellipticity gradients obtained from the PSPC data. We mention 
that the measured X-ray ellipticity gradients require steep mass profiles [p ~ r""*) for our SMD 
models; i.e. models where the mass is stratified on concentric spheroids of constant ellipticity. Of 
course, this does not exclude the possibility of a fiatter profile where the mass changes shape with 
radius. Models with small temperature gradients typical for rich clusters (Mushotzky 1995) affect 
the isothermal e^ass estimates by < 10%. 

We computed the masses of the X-ray-emitting gas and the total gravitating matter in 
spheres of radii (0.5, 1.0, 1.5)hgQ Mpc using the 95% statistical uncertainties of our gas and total 
mass models and the 90% uncertainties on the gas temperatures from Edge et al. (1990). The 
effects of the X-ray ellipticity on the gas masses is generally less than a few percent in agreement 
with the result of White et al. (1994). The ratio of gas mass to total gravitating mass of the 

3/2 1 

clusters increases with radius and has values Mgas/Mfot = (4% — ll%)hgQ within a l.ShgQ Mpc 
radius, in excellent agreement with the results of White & Fabian (1995) and thus consistent with 
the "Baryon Catastrophe" proposed by White et al. (1993). 

The shapes derived for the dark matter distributions using the isothermal models are 
essentially identical to those obtained for the total gravitating matter: eoM ~ ^tot + 0.02. The 
observed ellipticities of the galaxy isopleths (Carter & Metcalf 1980) are consistent with the 
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ellipticities of the dark matter. Moreover, the ellipticities of dark halos predicted by N-body 
simulations of a universe dominated by cold dark matter (Frenk et al. 1988; Efstathiou et al. 
1988; Katz & White 1993) agree with the results for the clusters in our sample. 

Our modeling procedure has generated tight constraints on the intrinsic ellipticities and 
density profiles of the dark matter from the X-ray data. To obtain the most general constraints 
a non-parametric estimation of the shape of the dark matter is warranted and should be feasible 
given the high S/N and spatial resolution of the observations along with the highly relaxed 
appearance of A401, A2029, and A2199 (see Buote & Tsai 1995b for a discussion of judging 
degree of relaxation). For these clusters a non-parametric treatment simply involves solving the 
hydrostatic equation for the potential in terms of the gas density and temperature which may be 
obtained from general spheroidal deprojection following Palmer (1994). 
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A. Prolate Ryden Potentials 

Ryden (1990) computes the gravitational potentials (up to quadrupole order) of oblate 
spheroids having ellipticity that varies with radius. Generalization of Ryden's result to prolate 
spheroids is a straightforward (albeit tedious) application of the multipole expansion of the 
three-dimensional gravitational potential (e.g., Binney & Tremaine 1987). Ryden considers the 
three-dimensional density distribution to be completely specified by the density, p{a), along the 
semi-major axis, a, and the eccentricity, e(a), of the isodensity surfaces. The isodensity surface for 
the prolate case is, 

r^cos^^ r^sin^^ 

+ a^l-e^) = ^' 

where r is the distance from the center, 9 is the conventional polar angle, and we have suppressed 
the dependence on a in e for notational convenience. Expressing the gradient in eccentricity as, 

„ a d(l — e^) , . , 

the mass of a thin shell of matter of width Sa and uniform density p is, 

5M = 4:k{1- e'^)pa'^5a[l + J:/3]. (A3) 
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The oblate potential is completely specified by the coefficients for the interior and exterior 
potential given by Ryden's equations (7) and (12). For the prolate case we obtain the following 
values, 



where, 
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Substituting the coefficients into equations (16) and (17) of Ryden gives the prolate gravitational 
potential up to quadrupole order. 
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Table 1: Cluster Observations 









O CV 1* /~\ 1 1 Tl /H 

JJdCKy,! U UliU. 


Oliistcr 


oec[uence ^ 


li/xposure (sj 


(^10 ^ cts s ^ Rrcmin ^) 


A /I ni 


rpoUU/oO 


1 400 






WpoUUloZ 


R7Q7 




A1656 


rp800005 


22183(20032) 


4.08 




rpSOOOOe 
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rp800009 


20691(19604) 






rp800013 


22427(20954) 




A2029 


rp800249 


12550 


5.50 


A2199 


wpl50083 


10563 


3.14 


A2256 


wplOOllO 
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2.50 




wp800162 
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wp800163 


10803(6690) 






wp800339 


4978(2437) 






wp800340 


9430(7119) 






wp800341 


10480(7469) 





Note. — Time-filtered exposure times are given in parentheses. 



Table 2: Pixel Scales 







Radial Profile 


■^max 


Ellipticity 


0.5-2.0 keV Flux 


Cluster 


z 


(arcsec) 


(hgo kpc) 


(Ko kpc) 


(arcsec) 


(hgo kpc) 


(10--^^ erg cm~^ s~^) 


A401 


0.0748 


40 


54.4 


1251 


20 


27.2 


29.6-29.7 


A1656 


0.0232 


105 


44.3 


487 


45 


19.0 


153 


A2029 


0.0768 


40 


55.7 


1448 


20 


27.9 


40.6-41.0 


A2199 


0.0299 


90 


48.9 


1321 


30 


16.3 


63.8-64.0 


A2256 


0.0581 


60 


63.4 


1458 


30 


31.7 


34.0-34.4 



Note. — Dmax is the edge of the outermost bin of the radial profiles. The fluxes are computed in a circle of radius 
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Table 3: (3 Models 



Cluster 


ax 
bf 


(Ko kpc) 
95% 


bf 


/? 

95% 




dof 


A401 
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173-185 


0.606 


0.597-0.617 


27 


20 


A1656 


267 


262-273 


0.785 


0.771-0.799 


136 


8 


A2029 


157 


141-174 


0.682 


0.664-0.700 


27 


16 


A2199 


85 


77-92 


0.653 


0.643-0.663 


56 


17 


A2256 


397 


367-427 


0.908 


0.870-0.947 


53 


9 



Note. — 



The 95% confidence limits reflect only statistical uncertainties on two interesting parameters. 
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0, 


,180-0, 


,210 


42.9 


39.6-46.1 




326 





.178 


0, 


,175-0, 


,200 


40.8 


37.7-44.0 




408* 





.164 


0, 


,150-0, 


,180 


36.4 


33.2-39.5 




489 





.149 


0, 


,135-0, 


,165 


36.8 


33.1-40.7 




571 





.146 


0, 


,125-0, 


,160 


34.3 


30.2-38.1 




652* 





.142 


0, 


,115-0, 


,155 


34.3 


29.4-38.8 


A2256: 


(ain = 317) 


















634 





.286 


0, 


,280-0, 


,290 


122.6 


122.1-123.1 




792 





.243 


0, 


,235-0, 


,250 


123.6 


122.8-124.3 




951* 





.211 


0, 


,205-0, 


,215 


123.2 


121.9-124.4 




1109 





.197 


0, 


,190-0, 


,205 


122.1 


120.9-123.3 




1268 





.185 


0, 


,170-0, 


,185 


120.7 


118.9-122.4 




1426* 





.164 


0, 


,130-0, 


,155 


119.3 


116.6-121.4 



Note. — 



and Qout are respectively the inner and outer semi-major axes of the aperture. 
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Table 5: 


Shapes of The 


Total Gravitating Matter 


Cluster 


Model Oblate £mass 


Prolate (mass 


A401 


C (ry = 2) 


0.52-0.61 


0.48-0.55 


A1656 


C(ry=l) 


0.40-0.41 


0.37-0.38 


A1656 


C (ry = 2) 


0.42-0.43 


0.39-0.40 


A2029 


H (n = 3) 


0.46-0.49 


0.42-0.44 


A2199 


H (n = 3) 


0.56-0.64 


0.49-0.55 


A2256 


C (ry = 2) 


0.36-0.38 


0.34-0.36 



Note. — These are the 95% confidence uncertainties on tmass for the isothermal mass models of §3.2.1. The models 
are C = "cored" power law and H = Hernquist (see §3.1). 



Table 6: Effects of Temperature Gradients on Mass Shapes 



T Model Oblate emass Prolate f; 



p = 


-0.27 


0.51-0.56 


0.47-0.50 


p = 


-0.65 


^0.6 


^0.6 


p = 


+0.07 


0.43-0.46 


0.40-0.42 


p = 





0.46-0.49 


0.42-0.44 



Note. — These are the 95% confidence uncertainties on Cmass for A2029 for different temperature gradients (§3.2.2) 
where p corresponds to the temperature index in equation [5]. All of the models have Hernquist (n = 3) densities 
(see §3.1) and the p = temperature index model is simply the isothermal case given in Table 5. 



Table 7: Gas Mass and Composite Mass (10^'* Mq) 



Cluster 


0.5/1 

Mgas 


8o' Mpc 

Mtot 


% 


l.Oh 

Mgas 


8o' Mpc 

Mtot 


% 


1.5/1 

Mgas 


8o' Mpc 

Mtot 


% 


A401 


0.276-0.288 


9.2-10.6 


3 


0.863-0.884 


14.1-25.0 


3-6 


1.54-1.64 


16.0-35.8 


4-10 


A1656 


0.258-0.259 


4.4-7.5 


3-6 














A2029 


0.287-0.308 


6.0-13.0 


2-5 


0.799-0.843 


11.4-22.5 


4-7 


1.38-1.53 


14.5-26.8 


5-11 


A2199 


0.144-0.147 


9.4-12.7 


1 


0.387-0.408 


12.2-17.0 


2-3 


0.657-0.711 


13.2-18.3 


4-5 


A2256 


0.232-0.256 


3.4-5.9 


4-8 


0.725-0.753 


9.4-16.8 


4-8 


1.16-1.28 


12.9-22.3 


5-10 



Note. — These masses are computed in spheres of radii (0.5, 1, 1.5)hgg Mpc. The "%" column denotes the percent 
ratio Mgas I Mtot. The gas mass uncertainties represent the 95% errors of the /3-model parameters and the 90% errors 
in the temperatures from Edge et al. (1990) as they affect determination of the X-ray flux and plasma emissivity. The 
composite mass uncertainties represent the 95% statistical uncertainties from the mass models and the 90% errors in 
the temperatures from Edge et al. (1990). 
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Table 8: Shapes of The Dark Matter 



Cluster 


Mo 


del 


% Mgas/Mtot 


Oblate Com 


Prolate £dm 




(-gal 


A401 


C (ry 


= 2) 


10 


0.53-0.62 


0.48-0.56 


0.625 


0.46-0.66 








20 


0.54-0.63 


0.49-0.56 






A1656 


C (ry 


= 2) 


5 


0.43-0.44 


0.40-0.41 


0.313 


0.51-0.71 








10 


0.43-0.44 


0.41-0.42 


0.625 


0.40-0.54 


A2029 


H (n 


= 3) 


10 


0.47-0.50 


0.42-0.44 


0.625 


0.41-0.61 








20 


0.48-0.51 


0.43-0.45 






A2199 


H (n 


= 3) 


5 


0.57-0.65 


0.50-0.56 


0.625 


0.47-0.67 








10 


0.58-0.66 


0.51-0.58 


1.25 


0.45-0.63 


A2256 


C (ry 


= 2) 


10 


0.37-0.39 


0.34-0.36 


0.469 


0.50-0.70 








20 


0.39-0.41 


0.34-0.36 







Note. — These are the 95% confidence uncertainties on toM for the isothermal mass models of Table 5 where now 
the models correspond to the dark matter for different ratios of gas mass to total mass, rgai = \/ab is the geometric 
mean radius (in h^g Mpc) of the aperture used to compute the galaxy isopleth ellipticities, Cgai, by Carter & Metcalfe 
(1980). 
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Fig. 1. — Contour plots of the reduced images for each cluster which have been smoothed with a 
circular gaussian {a = 15") for display purposes only. The contours are separated by a factor of 2 
in intensity and the displayed regions correspond to slightly larger than the region analyzed (i.e. 
out to radius Dmax)- The lowest contour level shown for each cluster corresponds to intensities (in 
10"'* counts cm-2 s'^) of 1.56 for A401, 40.8 for A1656, 1.77 for A2029, 5.05 for A2199, and 5.00 
for A2256. 

Fig. 2. — Shown are the background-subtracted radial profiles for each cluster. The top axes are 
labeled in hg^ kpc and the error bars are placed at the center of each radial bin and the boxes 
indicate bins omitted from analysis. 
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